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We investigate stability of two branches of Freund-Rubin compactification from thermodynamic 
and dynamical perspectives. Freund-Rubin compactification allows not only trivial solutions but 
also warped solutions describing warped product of external de Sitter space and internal deformed 
sphere. We study dynamical stability by analyzing linear perturbations around solutions in each 
. . . branch. Also we study thermodynamic stability based on de Sitter entropy. We show complete 

0^ ' agreement of thermodynamic and dynamical stabilities of this system. Finally, we interpret the 

, results in terms of effective energy density in the four-dimensional Einstein frame and discuss cos- 

' mological implications. 
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' ' I. INTRODUCTION AND SUMMARY 

De Sitter or quasi-de Sitter spacetimes describe the inflationary epoch of the universe at its early stages, and also 
_ the present universe which has entered the period of accelerated expansion. One of the most intriguing issues today 
is to realize such de Sitter vacua in fundamental high-energy physics. Particularly, higher-dimensional spacetimes 
^ are required by string theory, which is a promising candidate for the fundamental theory. In order to obtain the 
'— ' . effective four-dimensional theory in higher-dimensional spacetime we should usually compactify the extra dimensions 
■ and stabilize the compactified internal space. Thus we need to embed a four-dimensional de Sitter spacetime into 
^ , higher-dimensional spacetimes together with stabilization of the extra dimensions. 

' The Freund-Rubin compactification [l[ is a simple model with a stabilization mechanism that the extra dimensions 
QQ ' are dynamically compactified and stabilized by a flux of anti-symmetric tensor field or form field. In this model 
f — ■ with [p -\- q)-dimensional spacetime, a g-form flux field is introduced to stabilize the g-dimensional compact space. 
' Moreover, turning on a positive bulk cosmological constant allows an external de Sitter space and an internal manifold 
. with positive curvature (also, adding a dilaton field Q). Consequently, we obtain a (p + (7)-dimcnsional product 
' spacetime which consists of a p-dimensional de Sitter space dSp and a q-dimensional sphere S'', that is the Freund- 
0^ Rubin (FR) solution. 

It has been known that the Freund-Rubin solution has two classes of dynamical instabilities 0, Q ; one is attributed 
J> to homogeneous excitation {I ~ mode) of the internal space, corresponding to change of the radius of the extra 
r '~j ■ dimensions; and the other is inhomogeneous excitation with quadrupole moment (1 — 2 mode) and higher multi-pole 
rS ^ moments (/ > 3 modes), representing deformation of the extra dimensions. 

^ . The I — mode is the so-called volume modulus or radion, and becomes tachyonic when the Hubble expansion rate 
■ " " ' of the external de Sitter space is too large (in other words, when the flux density on the internal space is small). This 
fact implies that if the energy scale of the inflationary external spacetime is sufficiently larger than the compactification 
scale of the internal space, the volume modulus will be destabilized. In order to avoid the emergence of this instability, 
configurations with small Hubble expansion rate is preferable. 

While instability of the volume modulus exists already in g = 2, instabilities arising from deformation of extra 
dimensions emerge only if the number of extra dimensions is larger than or equal to four. In the unstable region, in 
which at least one of the / > 2 modes are tachyonic, the external spacetime has small Hubble expansion rate including 
the Minkowski spacetime. It means that the flux densities are very large in the unstable region of this type. It should 
be noted that for more than four extra dimensions, the two unstable regions overlap so that stable configurations for 
the FR solution no longer exist. 
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Little has been known about the non-perturbative properties of the instabihty from higher multi-pole modes;-'^ how 
it turns out after the onset of this instability and whether any stable configuration exists as a possible end-state in 
this model, and so on. In the previous work [0], one of us has shown that in the Freund-Rubin compactification there 
is a new branch of solutions other than the FR solutions. Those solutions are described as the warped product of an 
external de Sitter space and an internal deformed sphere. It has been found that the branch of the FR solutions and 
that of the warped solutions intersect at the point where the FR solution becomes marginally stable for the I = 2 
mode. Although we have seen existence of non-trivial solutions other than the FR solutions, their stability remains 
unanswered. 

In this paper we are particularly concerned with the close connection between dynamical stability and thermody- 
namic stability. The interesting relationship between dynamical and thermodynamic stability, which is well known 
as the correlated stability conjecture (or the Gubser-Mitra conjecture B, B ) , has been su gge sted and confirmed for 

some black objects (strings, branes and so on) by many authors I, [illiiriiS, M M MM [13 ■ (See e.g. [isl, [igf 

and references therein.) It is important to examine whether such connections really exist and whether they can be 
extended to systems other than black objects such as spacetimes with de Sitter horizons. 

In fact, for the FR solutions we can simply reinterpret the instability from the I — mode based on thermodynamic 
arguments ^] as follows (and also see ^]). For a fixed total flux on the internal space, the branch of FR solutions 
are divided into two sub-branches in terms of entropy defined by the total area of the de Sitter horizon: a sub-branch 
of solutions with higher entropy and the other with lower entropy. Both sub-branches terminate at one critical point, 
where the I = mode becomes massless and the FR solution is marginally stable. Moreover, the solutions on the 
lower-entropy sub-branch, which are thermodynamically unfavorable, are dynamically unstable since the I = mode 
is tachyonic. Thus thermodynamic instability exactly coincides with dynamical instability for the / = mode of two 
sub-branches within the FR branch. 

The aim of this paper is to examine the stability of the new branch of warped solutions from both dynamical 
and thermodynamic perspectives. This opens up new possibilities of the applicability of close connection between 
dynamical and thermodynamic stabilities. 

We examine the dynamical stability by analyzing perturbative stability of the system in a straightforward way. For 
simplicity, we restrict our considerations to the sector which behaves as scalar with respect to the external de Sitter 
space since unstable perturbations of the FR solutions are in this sector. 

In the case of four-dimensional external spacetime and four-dimensional internal space, we numerically obtain the 
Kaluza-Klein (KK) mass spectrum and show that the warped solutions are stable in the low Hubble regime, while 
the FR solutions are unstable due to the I = 2 mode in the same regime within numerical accuracy. 

In order to reveal the thermodynamic property we derive the first law of de Sitter thermodynamics for Freund- 
Rubin compactifications. Each branch of the FR solutions and the warped solutions obeys the first law in terms of 
entropy S and total fiux $: 

as ^ -Jhiz^d<^^ (1) 

4(p - l)hP ^ ' 

where two parameters h and h characterize the flux density on the internal space and the Hubble expansion rate of 
the external de Sitter space, respectively. This fact means that for a fixed total flux, the branch with higher entropy 
should be thermodynamically favored. Comparison between the entropy of the FR branch and that of the warped 
branch for a given total flux tells us which branch is thermodynamically favored. The result is that the warped branch 
is entropically favored in the low Hubble regime while the FR branch is favored in the high Hubble regime. 
The above results are briefly summarized as follows: 

• For small Hubble expansion rate, the warped branch is thermodynamically favored and dynamically stable. 

• For large Hubble expansion rate, the FR branch is thermodynamically favored and dynamically stable. 

Thus, as we have expected, we see complete agreement of thermodynamic and dynamical stabilities for two branches 
of FR compactifications. This provides yet another example showing close connections between thermodynamic an 
dynamical properties of systems with horizons. 

The rest of this paper is organized as follows. In Sec. [Til we review general Freund-Rubin compactifications and 
show the Freund-Rubin solutions describing dSp x S'^ and the warped solutions describing a warped product of an 
external de Sitter space and an internal deformed sphere. In Sec. Illll we investigate dynamical stability of the warped 
solution by considering perturbations around the background solution. In Sec. IIVI we derive the first law of de Sitter 
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thermodynamics and discuss thermodynamic stabihty for the FR branch and the warped branch. In Sec. |V] we 
interpret the above results in terms of effective energy density in the four-dimensional Einstein frame and discuss 
cosmological implications. 



II. FREUND-RUBIN COMPACTIFICATION 



In this section we review general Freund-Rubin flux compactifications, including a bulk cosmological constant. We 
consider the {p + (7)-dimensional action 

where A is a (p + (7)-dimensional bulk cosmological constant and is a q-form field strength for stabilizing the 
q-dimensional internal manifold. (We use units in which G = I imless otherwise noted.) The Einstein equation and 
the Maxwell equation lead to 

GmN = _ ^y PMLi-'-Lg-iFN^^'"^"'^ - -^^P(q)9MN - ^QMN (3) 

and 

y^^pMN,-N,^, = 0, (4) 

where the 5- form field satisfies the Bianchi identity V[Mi-FM2 - Mg+i] = 0- 



A. Freund— Rubin solution 

These equations have well-known solutions originally found by Freund and Rubin T|. The metric and the g-form 
field strength in the Freund-Rubin solutions are given by 

ds^ = -dt^ + e^'''dxl_i + p^dn^, (5) 

and 

= k^^...^^, (6) 

where e^^i-..^ is the volume element of the g-sphere with a radius p. These solutions describe the direct product of 
a p-dimensional external de Sitter space with Hubble expansion rate h and a g-sphere with radius p. The internal 
g-sphere is supported by the g-form flux with a flux density b. The Einstein equation and Maxwell equation yield 
relations among b, h and p: 

{p-l)ip + q-2)h^ + iq-l)b' = 2A, (7) 

{q-lfp-^ + ip-lfh^ = 2A. (8) 
Evidently, these relations allow a one-parameter family of solutions. 



B. Warped solution 



There is another one parameter family of non-trivial solutions. Their geometry is a warped product of de Sitter 
space and a deformed sphere. Such warped solutions are described by the following ansatz for the metric 

ds^ = e2^M[-dt2 + e2'^*df2_J + e-F&'^('-)[dr2 + a^{r)dnl_,], (9) 

and the g-form flux 

= 6e-^'^'^a«-Mr A df),_i, (10) 
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FIG. 1: Two branches of solutions in the {h^ ,h?) plane. The blue solid line represents the branch of Freund-Rubin solutions. 
The red bold points correspond to values calculated numerically for which warped solutions have been found. Two branches 
intersect at one point (b^ ,h^) — Note that we have set A = 1. 




where 6 is a constant. Note that (|10p automatically satisfies the Maxwell equation and the Bianchi identity. Then, 
from the Einstein equation we have the two equations, 



9 , p{p-l] 



a q — 2 



-e 



q-2 



2A _ 2p ^ 
-e 9-2'' 



and the constraint equation 



+ 



p + q-2^ 



q-2 

2p(.?-i) J. 2 A 



(9-2) 



a — 1 



p + q-2 



e-^^-iq-l)%', 
a 



{q-l){q-2) 



1 



p{p + q-'2),,2 , Pb-1),2 -2(£±s^0 , -Ma 
2{q-2) +2"^ 



__2p_ 

Ae 9-2^ 



(11) 



(12) 



We are interested in the case where the internal space is compact. Thus, we consider an interval r_ <r<rj^, where 
a{r) vanishes at the endpoints, a{r±) = 0, and is positive between them. Then the regularity requires the following 
boundary conditions: 



|a'(r-±)| = l, 0'(r±) = O 



(13) 



at the endpoints. These conditions ensure that the internal space with spherical topology is regular at the north and 
south poles. 

A one-parameter family of solutions to the above equations and boundary conditions for p = A and q — A was 
found numerically by one of the authors [H] ■ This warped branch of solutions emanates from the marginally stable 
solution in the branch of Freund-Rubin solutions as shown in Fig [1] In Fig [1] two lines representing two branches 
intersect at one point. At the intersection, the solution in the warped branch is no more warped and reduces to a FR 
solution. For &^ smaller than the value at the intersection, the internal space is prolate. On the other hand, for 
larger than the value at the intersection, the internal space is oblate. Therefore, while the numerically obtained value 
of at the intersection includes numerical errors, the statement that the warped solution reduces to the FR solution 
at the intersection is exact: there must be a boundary between the oblate and prolate cases; and at the boundary, 
the internal space is spherical and the warp factor is trivial. 

In the following sections, we shall investigate stability of the FR and warped branches. 



III. DYNAMICAL STABILITY 



In this section we investigate dynamical stability of FR solutions and warped solutions by considering linear per- 
turbations around them. We will concentrate on scalar-type perturbation with respect to the external de Sitter space 
since in the case of the FR solution instability arises from perturbations of this type. 
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A. Background 

In the previous section we have aheady shown the background metric and form field. In this subsection we rewrite 
them in a form which is convenient for the analysis of perturbation equations. 
We suppose that the {p + (7)-dimensional metric is given by 



A\y)g^.,{x)dx''dx'' + dy^ + B^{y)j,j{z)dz'dz\ 



(14) 



where g^i, is the metric of p-dimensional Lorentzian Einstein space and jij is the metric of (q— l)-dimensional Euclidean 
Einstein space. Note that the Riemann tensors with respect to g^i, and 7.^ respectively satisfy '^^'^R^u = K{p— l)g^i, 
and ^'^^Rij — k{q — 2)7^ . The q-form field strength is 



q-l 



(15) 



which satisfies the Maxwell equation and the Bianchi identity automatically. 

Then, non- vanishing components of the (j) + (7)-dimensional Einstein tensor Gmn become 



2 + + 2 W 



k , B" A'B' 

- + ('?-i)^ + b-i)('z-i)^ 



A' 



p{p-l)A''~K , {q-l)iq-2)B''-k , A'B 

Gyy^ ^ + ^ E^+Piq-i) 



A^ 
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AB ' 



Gij — 



p(v-l)A'^-K A" iq - 2){q - i) B'^ - k , „ B" , ^A'B' 

2 A^ 2 B^ +(-^-^)^+^('^-^)af 



(16) 

(here in this section the prime denotes the derivative with respect to y) and the energy-momentum tensor of the 
q-form fiux field is given by 



T — — 



9 

2^2(p-i)5m^' 



Tyy - 2A^p ' 



(17) 



and the other components vanish. 

Finally, we note that if we set A = e*^, B = e~'^'^a, dy = e~'^'^dr, fc = 1 and K ~ in the above equations, 
then all equations for the background ansatz in Sec. are reproduced. 



B. Perturbation 



We now leave the subject of background and turn our attention to linear perturbations around the background. 
We suppose that the p-dimensional external spacetime and the (g — l)-dimensional internal space are the Einstein 
manifolds with the metric and 7^ , respectively. In this case we can decompose tensors on the (p + g)-dimensional 
spacetime into scalar- type, vector- type and tensor- type components with respect to 5^,^ and jij. Hence, in our analysis 
we decompose perturbations of the metric and the form field into different types and obtain decoupled perturbation 
equations in each sector. 

As we already mentioned before, we will concentrate on scalar-type perturbations. Especially, we suppose that the 
internal manifold is topologically a sphere and has SO((7)-isometry, namely jij{z)dz^dz^ — dfl^_^ which is the metric 
of unit round {q — l)-sphere. Then what we are interested in is perturbations which are scalar- type quantities with 
respect to not only p-dimensional de Sitter symmetry of the background external space but also SO{q) symmetry of 
the background internal space. 

In this paper, for simplicity we assume that the perturbations preserve the SO(q) symmetry of the background 
internal space. By choosing an appropriate gauge (see Appendix [X]) . we can write the perturbed metric and field 
strength as follows: 

ds^ = (1 + BY) A2(y)g^,dx^da;'' + [1 + (B - n)Y] dy^ + 



n 



^] 

q-l q-l 



B Y 



(18) 
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and 



(9) 



59-1 

&^^dy A dnq^i + d (^Y) A dOg_i 



(19) 



where Y(a;) is scalar harmonics on the de Sitter space with the Hubble expansion rate h and variables 11, O and 
depend on only y-coordinate due to the SO(Q')-symmetry. 

From the linearized Einstein equation and Maxwell equation, we obtain a set of two perturbation equations for two 
variables 11 and Q: 



{p + q- 2)n" + {q- 2)fl"+{p + q - 2) 



n' + {q- 2) 



f2' 



[(p + (7-2)n-gn]=0, 



0" + 



(3p-2)^+3(g-l)|^ 



+ 



, 2(g-2) 

2(p + g~2)(g-2) 
m2 + 2/i2(p-1)2 2(;(g-2) 



-4A 



(20) 



n 



^2 



+ 
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-4A 



= 0, 



where /x^ is the KK mass squared which is defined by V^Y(a;) = \x^y{x). Here denotes the covariant derivative 
with respect to g^^v Moreover we have an algebraic equation for y>: 



LP : 



26 



(21) 



Boundary conditions are specified by the regularity at the poles of the internal space of spherical topology, which are 
characterized by -B = 0. They are given by 



{p + q-2)Il-qQ = fl' = at B = 0. 



(22) 



Thus the perturbation equations are reduced to eigenvalue problems with eigenvalue fi^. If the spectrum of fi"^ is 
non-negative, we can conclude that the background spacetime is dynamically stable. 

Before we discuss stability of the warped solution, let us recall stability of the FR solution. For the FR solution, 
we set A = 1 and B = p sin ^ . Then the perturbation equations reduce to 



{p + q- 2)H" + {q- 2)il"-{p + q - 2){q - l)^n' + {q - 2){q - 



+ 



B' 



n" + 2,{q-l)—n'- 
B 



+ 



2 , 2(g-2) 
^ + 52 
2(p + g-2)(g-2) 

B2 



[(p + g-2)n-gfi]=0, 



4A 



(23) 



H 



m' + 2/i2(p-1)2 + M|_^_4A 



f2 = 0. 



Eliminating O from these equations we obtain a single forth order differential equation for 11, 



£)2 . £)2n + 2 



p + q-2 



H = 0, 



(24) 



where denotes Laplacian on 5^ with a radius p. We expand H in terms of scalar hamonics Y[y) on 5^, and then 
the mass eigenvalues p? are given by 



,2 , («-l)(^'-2)^2_(p_i^^2 



p + q-2 



± 



(g-l)(p-2) 
p + q-2 



n 2 



b^-{p- l)/l2 



p + g- 2 



(25) 
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where D'^Y{y) = -\Y{y) with A = l{l + q-l)p-'^ It is clear from this expression that the scalar perturbations 

for each multi-pole moment I generically have two independent modes corresponding to the mass eigenvalues /i^ and 
/i'i. However, for / = and / = 1, only one of them is physical and the other is a gauge mode. For / > 2, we denote 
the modes with the mass squared /i^ and /i^ as / — 2+, 3+, • • • and I = 2_, 3_, • • • , respectively. In the I = 1 case 
the physical mode has the mass eigenvalue = 1). On the other hand, in the I — case we have a physical 

mode with ^^{1 = 0) = ~ 2(p - l)h^, corresponding to nl for ^'^^^^^b^ > {p ~ l)h^ and ^il for 

Thii 
where 



This mass spectrum leads to the following result for dynamical stability of the FR solutions; when > h? 



c(i=0)' 



_ 2A(p-2) 

/ic(;=o)- (p_l)2(p + ^_2)' ^^^^ 

the I = mode is tachyonic and the FR solution is unstable arising from homogeneous excitation of the sphere. In 
addition, for g > 4, when becomes smaller than the critical value /ic(/=2) by 

,2 _ 2A[(p~l)g^-(3p-l)g + 2] 

q{q-i){j>-mp + q-2) ' ^''^ 

the mass squared /i^ is negative for I = 2 and the FR solution is unstable arising from inhomogeneous excitations. As 
we have explained above, we call this mode 1 — 2^. Some modes with I > 3 can be unstable when is even smaller. 
The mass spectrum for the scalar perturbations of the Freund-Rubin solution is shown as blue dashed lines in Fig [21 

Now, let us consider the warped solutions. Especially we would like to examine dynamical stability of the warped 
solutions in the small Hubble regime in which the FR solutions suffer from instability of the inhomogeneous excitations. 
For the warped solution we numerically solve the eigenvalue problem for a set of differential equations (j20p with the 
boundary condition (|22p in order to find the mass spectrum /i^. In Fig [21 we present the numerically obtained mass 
spectrum of the warped branch for p — A and g = 4. The red solid lines represent for the warped branch and the 
blue dashed lines for the FR branch. In the left panel of Fig [21 we focus our attention to the / = 2_ mode since it 
is the first unstable mode for inhomogeneous perturbations on the FR branch. Note that "Z = 2" means quadrupole 
moment with respect to the SO(g-|- l)-symmetry of the internal space and that, rigorously speaking, this terminology 
is valid only for the FR branch. However, as explained in the end of Sec. [Hi there is a critical value of 6^ at which 
the solution in the warped branch reduces to an unwarped, FR solution. Since h? is determined by 6^ in each branch, 
this implies that there is a critical value of h? at which the solution in the warped branch and that in the FR branch 
represent the same solution. Therefore, it makes perfect sense to define the "Z — 2 mode" for the warped branch 
as the mode which approaches the I = 2 mode of the FR branch as approaches the critical value. We can 

define I = 3±, 4±, • • • modes for the warped branch in a similar way. As one can see in the left panel of Fig [21 for 
h? < /ic(;=2) ~ A/18, /i^ for the / = 2_ mode of the warped branch is positive, while the I = 2_ mode of the FR 
branch becomes tachyonic. Thus the warped branch is stable configuration in the low Hubble region, in which the 
FR branch is unstable. Actually, within numerical accuracy one can see that the red solid line and the blue dashed 
line for I — 2^ intersect at = 0. This implies that the critical value of at which the warped solution reduces to 
the FR solution agrees with (|?7|) . 

The mass squared for some other modes as well as the I — 2_ mode is shown in the right panel of Fig [21 As seen 
from the mass spectrum, the warped branch has no unstable mode in the low Hubble region. In addition we notice 
that /i^ of the warped branch (the red solid line) is larger than that of the FR branch (the blue dashed line) for 
smaller than /i^(;^2) (^^^ ^^^^ hand side). This means that deformation of the internal space and warping tend 
to stabilize the shape modulus of the internal space in the low Hubble region. In other words, the tachyonic shape 
modulus is stabilized by the condensation of the modulus itself. 



IV. THERMODYNAMIC STABILITY 



In the previous section we have investigated dynamical stability of two branches of Freund-Rubin compactification. 
In this section we shall investigate thermodynamic stability of the same system and compare the results with dynamical 
stability. 
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FIG. 2: The mass spectrum for scalar perturbations: the Z = 2_ mode (left) and some modes with Z = 0, 1, 2±, 3± (right). Red 
solid lines and blue dashed lines indicate the warped branch and the FR branch, respectively. The green vertical line indicates 
the critical value — /i^^j^j) = A/18 at which two branches merge. In the low Hubble region where the FR branch is unstable, 
the warped branch is stable. 



A. Thermodynamic relations 



In Freund-Rubin compactifications we can define various physical quantities characterizing thermodynamic prop- 
erties of the system. One of the most important among them is the de Sitter entropy S, which is defined by one 
quarter of the total area A of de Sitter horizon. For the metric ([9]) it is given by 



_ ,4 _ rip_2^g-i 
" T " ihP-^ 

Also, we can define the total flux of the (/-form field (fTU]) as 



are ^a^ . 



(28) 



dre <j-2 Vq; 



9-1 



(29) 



which is a conserved quantity for this system. 

Before examining thermodynamic stability, let us derive the first law of de Sitter thermodynamics for Freund-Rubin 
compactifications. 

For this purpose it is convenient to consider variations of the Euclidean action for the system since the on-shcU 
Euclidean action is directly related to the de Sitter entropy as we shall see below. Assuming SO(p + 1) x S0((7) 
isometry, we can take the metric ansatz as 



(30) 



where dflp and dfl^_i denote the metrics of the unit round p- and {q — l)-sphere, respectively. The g-form field 
strength is given by 



The Euclidean action is given by 



'Euclid 



^/d^+^.EV^(i?-2A- 



p2 



(31) 



(32) 



The Ricci scalar and the field strength are 



2p 

R =e~'' 



,'2 



('?-l)('?-2)- 



+ 1 p{p + q-2)^, 



2p 



,9-1 



q-2 

p J./ 



(33) 



a q — 2' 



9 



and 



2pq 



Hence we have the foUowing expression for the EucUdean action 



^Euciid[a, 0, "0] 



WnhP 



2p 



dr 



{q-l){q-2)- 



1 



pi:? 



(34) 



b' +p{p-l)h^e ^ 



e "-^ ,2 

02(9-1) 



(35) 



,9-1 



where the boundary terms have vanished since the boundary conditions require a(r±) — 0, |a'(r±)| = 1 and (j)'{r±) — 
at the boundaries r = r±. Using the equations of motion, we evaluate the Euchdean action /Euclid on shell and it 
turns out that 



^Euclid — ~S, 

as explained in appendix |B] Also, the equation of motion for the form field is given by 



e 

a9-i 



= 0, 



and it can be easily integrated as 



e 



a9- 



— V' = b, 



where b is an integration constant. The total flux <f> is rewritten as 



2p(g-l) . 



(36) 



(37) 



(38) 



(39) 



We are now ready to derive the first law of de Sitter thermodynamics in our setup. We consider the first variation 
of the action /Euclid [a, 0, "0] with respect to a, (f) and tp- Suppose both {a,(p,ip} and {a + 5a, (j) + Scjy.ip + ^V"} are 
different sets of solutions satisfying the equations of motion, the first variation of the action (5/Euciid is given by 

2M9-i) ^ r+ 



5h 



Euclid 



ai- 



dr (EOM for a, (f) and V') • 



(40) 



The integrand in the last term will vanish because of the equations of motion for a, (f> and ip. As a result, only the 
boundary term contributes to the first variation of the action. By using ((38|) and ((39|) the first law can be derived as 

rip_2^ 



dS = -- 



-d$. 



(41) 



A{p- l)hP 

This implies that the entropy is described by a function of the total fiux. Hence the entropy S' is a thermodynamic 
potential with respect to the total flux <I> as a natural thermodynamic variable. We shall call a sequence of solutions 
which satisfles the first law a "branch" of solutions. 

So far, we have assumed that the bulk cosmological constant A is not a dynamical variable but a given constant. 
However it is probable that A is a dynamical variable induced by dynamical fields such as a scalar field. In this case 
we consider A as an additional thermodynamic variable so that the de Sitter entropy is now a function of <I> and A, 
S'(<I>, A). The first law (|4ip is easily generalized to include variation of A as follows. First, dimensional analysis leads 
to the following scaling relation 

5'(A-(9-i)/2$^^A) = A-(P+«-2)/25'($^A). 
Second, taking derivative with respect to A and setting A = 1, we obtain 

dS p+q~2^ 



q - I ^dS , , 

- $ h A — 

2 dA 

This is the first law of de Sitter thermodynamics for S{^, A): 



(42) 
(43) 



dS = - 



4(p - l)hP 



d$ - 



-S- 



- 1 ilp_2& 

'2~A{p-l)hP 



dA 



(44) 
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B. Stability 

In this subsection we discuss thermodynamic stabihty of the FR branch and the warped branch. As we have seen, 
the entropy S is the thermodynamic potential when we choose the total flux $ as a natural variable. Therefore, the 
second law of thermodynamics states that, for a fixed value of the total flux, a configuration with larger entropy is 
thermodynamically favored. 

To begin with, let us examine thermodynamic property of the FR branch. We shall see that the FR branch has 
two sub-branches and one of them is thermodynamically preferred than the other. In the FR branch, the entropy S 
and the total flux $ are given by 

^=%#l^> 'J> = ^V (45) 

where we have used Eqs. ([28]) and ([39]) with a{r) — psin ^ and (j){r) = 0. Note that we can explicitly check that these 
quantities satisfy the first law (|4ip . We find that the entropy is written as a function of the total flux and splits the FR 
branch into two sub-branches: a lower-entropy sub-branch and a higher-entropy sub- branch, as shown in Fig [31 (For 
example, see 20].) Therefore, for a given total flux the higher-entropy sub-branch is preferred than the lower-entropy 
sub-branch within the FR branch. 

The critical point dividing the FR branch into two sub-branches is determined as follows. As seen before, the FR 
branch satisfies the first law. However, the entropy 5'($) is a double-valued function of the total flux <I>. Nonetheless, 
the FR solutions can be described as an one-parameter family of solutions, for example, in terms of the Hubble 
expansion rate h of the external de Sitter space. Actually, the entropy is a single-valued function of h. These facts 
mean that a map from $ to ft, becomes singular at the critical value. Hence we can obtain the critical point by solving 
d<f>/d/i = dS/dh — 0, which yields the critical value 

2 _ 2A(p^2) _ (p-l)(p + g-2) 

'*^('=o)~ (p-l)2(p + q-2)' ^-('=0)- (p-2)(g-l) ''<'='y ^^^> 

By comparing with (|26p , it is easy to see that this thermodynamical critical point agrees with the threshold at which 
the I = mode becomes massless on the FR branch. In addition, the lower-entropy sub-branch is dynamically 
unstable against homogeneous {I = 0) excitation of the internal space. Therefore, we see complete agreement between 
thermodynamic and dynamical stability of the two sub-branches. 

In the previous paragraphs we have compared entropies of two sub-branches within the FR branch. We now compare 
entropies of the FR branch and the warped branch. For simplicity we shall consider the case with p — -i and g = 4 as 
an explicit example. 

We denote the entropy of the FR branch and that of the warped branch as Sfr and Sw, respectively. A difference 
between 5'w and Sfr for various values of the total flux is shown in Fig. [H It turns out that for $ < 32V37r^A~^/^, 
or equivalently, for 

< 4 (47) 

the warped branch has larger entropy and thus is thermodynamically stable. On the other hand, for $ > 32V37r^A~^/^, 
or equivalently, for > the FR branch has larger entropy and is thermodynamically stable. At the critical point 
where two branches merge the solution is marginally stable. 

It is worth noting that the thermodynamic stability investigated here agrees with the dynamical stability examined 
in Sec mil The FR branch has dynamical instability arising from inhomogeneous (/ = 2) excitation when < /ic(;=2)' 
where ft^(;^2) given by ([27]) . On the other hand, the warped branch is dynamically unstable for > /i^(;^2)- 
For p = 4 and q = 4, the critical value ft.^(;^2) f'-"' ^^'^ dynamical stability agrees with the critical value for the 
thermodynamic stability given in the right hand side of (j47p . Therefore, we again see complete agreement between 
thermodynamic and dynamical stabilities. 



V. COSMOLOGICAL IMPLICATIONS 



In previous sections we have investigated stability of two branches of flux compactification from thermodynamic and 
dynamical perspectives. One branch has higher symmetry SO(g-|- 1), where q is the number of extra dimensions, and 
corresponds to unwarped, Freund-Rubin solutions [l|. The other branch has lower symmetry SO{q) and corresponds to 
warped solutions found recently by one of the authors f6j . By fixing or scaling out the higher-dimensional cosmological 
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FIG. 3: The entropy of the FR branch and the warped branch as functions of the total flux $. Difference between entropies of 
the two branches is so smaU that two hnes are indistinguishable in this figure. See Fig. |4] for the difference. 



Sw~Sfr 
\ 



530 540 ' 5^"*-~_560 570 



FIG. 4: The difference between the entropy of the warped branch Sw and that of the FR branch Sfr- 



constant, each branch is parameterized by one parameter, either total flux $ of an antisymmetric fleld or the Hubble 
expansion rate h of the 4-dimensional de Sitter metric. We have seen that the unwarped branch is dynamically stable 
for h larger than a critical value /ic(i=2) but unstable for smaller values. On the other hand, the warped branch 
is dynamically unstable for h > ft.c(i=2) a-nd stable for h < /ic(;=2)- To investigate thermodynamic perspective, we 
have defined the total de Sitter entropy S as the 4-dimensional de Sitter entropy integrated over extra dimensions. 
We have shown that the dynamically stable branch, i.e. the unwarped (or warped) branch for h > /ic(i=2) (or 
h < /ic(i=2)j respectively), always has larger total de Sitter entropy than the dynamically unstable branch. Therefore, 
thermodynamic stability agrees with dynamical stability. 

In this section we consider cosmological implications of these results. For this purpose we shall first define the 
4-dimensional Einstein frame. For the (4 -I- g)-dimensional metric of the form 



GMNdX'^^dX 



(48) 



the higher-dimensional Einstein-Hilbert action is 



{Mi+q) 



2+q 



Ml 



A'^Xyf^n'^R[g\ 



(49) 



where M^+q and M4 are (4 -I- q)- and 4-dimcnsional Planck scales (in this section we have temporarily restored the 
Planck scales), and 



Ml 



(50) 
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Since ft in general depends on the 4-dimensional coordinates x^^ , the resulting 4-dimensional effective theory describing 

( F1 

gfil, is not Einstein but a scalar-tensor theory, ft is convenient to define the 4-dimensional Einstein frame gl^j by 

giV = ^'9^.., (51) 
in terms of which /4+g now includes the 4-dimensional Einstcin-Hilbert term: 

/4+, = ^ / d^xv/^i?[5(^^] + • • • . (52) 



Suppose that the 4-dimcnsional metric g^^i, represents a de Sitter spacetime with the Hubble expansion rate h. The 

abblc expansion rate 

He = n^^h, (53) 



corresponding Einstein frame metric glijy has the Hubble expansion rate 



Se = SnHdih-/ = 8^'(^4+,)^+^ , d'^y^qA'. (54) 



and the de Sitter entropy 

Actually, this agrees with the total de Sitter entropy ((28l) : 

Se = S. (55) 

For cosmological considerations energy density is more convenient than de Sitter entropy since the former can easily 
be extended to a general FRW universe. The effective energy density in the Einstein frame is 

PE = ZMlh^ = i I d''2/^/^^^') • (56) 



This is related to the total de Sitter entropy as 



Pe Stt 



2 



3Mi S 



(57) 



Therefore, the results of previous sections are restated in terms of the effective energy density pE in the Einstein 
frame as follows: dynamically stable branch, i.e. the unwarped (warped) branch for h > /ic(;=2) (^ < ^c(/=2))i always 
has lower pE than the dynamically unstable branch. This strongly suggests that a solution in the dynamically unstable 
branch should evolve to a solution in the dynamically stable branch. The latter solution is uniquely specified by the 
former since the flux <i> conserves. By this evolution, pE decreases. 

Moreover, the results of the previous sections suggest a new type of phase transition. Suppose that A is not the 
genuine constant but has a contribution from dynamical fields such as a scalar field. In this case A is expected to 
decrease while the flux <I> stays constant. (See Fig. [51) If we start with h > /ic(i=2) then the unwarped branch has 
lower energy density and is stable. Thus a solution in the unwarped branch should be realized initially. However, as 
A decreases, h also decreases and can reach the critical value /ic(;=2)- At that point, the stable and unstable branches 
merge. After that, for h < /ic(i=2)j solutions in the warped branch should be realized since this branch has lower 
energy density and is stable. This phase transition should be second-order since only one of the two branches is stable 
at a given time. 

Second-order phase transitions play important roles in cosmology. For example, the end of hybrid inflation is due to 
second-order phase transition. Thus, this kind of phase transition may provide a new way of realizing hybrid inflation 
in higher-dimensional theories. This possibility will be investigated in future publications. 



VI. DISCUSSION 



In this paper we have investigated stability of two branches of Freund-Rubin compactification from two perspectives; 
one is thermodynamic stability based on de Sitter entropy, and the other is dynamical stability of linear perturbations 
around the background solutions. 

We have analyzed linear perturbations around the warped solutions in order to examine dynamical stability of the 
solutions. The warped solutions are stable if the Hubble expansion rate of the external de Sitter spacetime is low 
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FIG. 5: The difference of the effective energy density between the warped branch and the FR branch with a fixed total flux. 



enough. In the same regime of the Hubble expansion rate, the Preund-Rubin solutions have instability arising from 

the I = 2 mode. It follows from what has been said thus far that deformation of the internal space and warping 
will stabilize unstable configurations. This may be considered as spontaneous breaking of the symmetry of internal 
space in the sense that less symmetric configurations are dynamically chosen. Actually, for the reason explained 
below, this phenomenon is natural from gravitational viewpoints. In cosmology it is well known that configurations 
with high matter density suffer from gravitational instability due to long-wavelength modes, namely Jeans instability. 
Jeans instability develops inhomogeneities in the universe and results in structure formation. In the case of FR 
compactification, Einstein equation tells that the external space with small Hubble expansion rate corresponds to the 
internal space with large flux density. Therefore, it is natural to expect that when the Hubble expansion rate is low 
enough, the flux distribution may become inhomogcncous due to analogue of the Jeans instability. 

In order to analyze thermodynamic properties, we have first derived the first law of de Sitter thermodynamics in 
terms of entropy and total flux. Sequence of solutions belonging to a branch obey this first law when those parameters 
characterizing the solutions change. Since the entropy is a natural thermodynamic potential for the total flux, for a 
given total flux, configuration with higher entropy should be favored thermodynamically. For p = 4 and q = 4 we 
have compared the entropy of the FR branch with that of the warped branch. There is a critical value of the Hubble 
expansion rate at which two branches merge and it is thus obvious that the entropies of the two branches agree at 
the critical value. For smaller Hubble expansion rate, the entropy of the warped branch is larger than that of the FR 
branch and thus the warped branch is thermodynamically favored. On the other hand, for larger Hubble expansion 
rate, the FR branch is thermodynamically favored. We found complete agreement of thermodynamic stability and 
dynamical stability. 

It is intriguing to see that correlation between thermodynamic stability and dynamical stability exists for Freund- 
Rubin flux compactiflcation. For a certain class of black objects the existence of such correlation has been known, that 
is the so-called Gubser-Mitra conjecture. This conjecture has been explicitly checked to hold for various black strings 
and branes. From what has been discussed above, it is probably natural to expect that the concept of correlated 
stability can be extended to a wider class of gravitating systems. 
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APPENDIX A: GAUGE 

In the background spacetime given by 

gMNdx^dx"" = A^{y)g^,dx^'dx'' + dy^ + B^{y)y,jdz'dz' , (Al) 

where gfj,^ and 7ij are respectively the metric of p-dimensional de Sitter space and {q — l)-dimensional round sphere, 
we consider scalar-type perturbations, 

Sgij = (oiDj - /j(LL) ^ B\jh^^\ (A2) 

5giy=Dih^^\ 

^9yy — ^yyi 
Sgtii = V^A??- 

Note that and Di denote covariant derivatives associated with g^i, and 7,j , respectively. 
These components transform as 

//(LL) ^ ^(LL) _ ^-2^(L)^ ^(Y) ^ ^(Y) _ ^-22^2^(L) _ ^^^^^ ^f(L) ^ ^(L) _ _ ^2(^-2^(L)y^ 

p A 

(7 — 1 -D 
hyy-hyy-2^y', ^ ^ ^ _ ^ (D _ ^ W , 

By setting 

^(L) ^ ^2^(LL)^ ^ ^(L) _ ^2(^(LL)y^ (^4) 

and assuming that the perturbations docs not depend on ^;*-coordinates, we can simplify the form of metric pertur- 
bations so that non-vanishing components are 

5gu.u = A^g^yH'^^\ 6gij = B'^-fijhP^\ Sgyy = hyy, (A5) 

and the other components vanish. In addition, using parts of the linearized Einstein equations for scalar-type pertur- 
bations we have an algebraic relation 

ip-2)H'^^'> + hyy + {q-l)h^''^ = 0. (A6) 

Finally, we expand the perturbations in harmonics Y(a;) on the p-dimensional de Sitter space 5^^, and then we 
define two variables n(y) and f2(y) as 



= UY{x), hyy = (n - n)Y{x), /I'Y) = ^^-^ - 



Y{x). (A7) 



APPENDIX B: EUCLIDEAN ACTION AND ENTROPY 

The Euclidean action is given by 



/Euciid[a,<^,^] = - J dr 



Mini), 



2p , e 9-2 ^ 9 
a2(9-i) 



a 



a2 q-2 

9-1 (Bl) 



Clp^q—l 

IG-rrhP 



J drjC{a,a',^,<p',i>'). 
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Here, since this action is invariant under the following transformation: 
where A is an arbitrary constant, we have an identity 



p dC dC 
-.a+ — 



„ 2(p + g-2) , 

pC + 2p{p - l)h^e F^"*a«"i = 0, 



(B2) 



(B3) 



q-2da' 

provided that air), <f>{r) and ?/;(r) satisfy the equations of motion. Hence the on-shell Euclidean action can be written 



Euclid — 



^p^g— 1 



dr(p- l)e ^^-^a-J-i = -S, 



(B4) 



where the boundary terms vanish because of the boundary conditions. It is easy to show by using (p— l)Qp — 27rf2p_2 
that S agrees with the de Sitter entropy defined in ([28]) . 
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